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■ Abstract. Some OB stars show variable non-thermal radio emission. The non-thermal emission 

' is due to synchrotron radiation that is emitted by electrons accelerated to high energies. The 

, electron acceleration occurs at strong shocks created by the collision of radiatively-driven stellar 

' winds in binary systems. Here we present results of our modelling of two colliding wind systems: 

I 1 ; Cyg 0B2 No. 8A and Cyg 0B2 No. 9. 

' Keywords, plasmas, radiation mechanisms: nonthermal, methods: numerical, binaries: spec- 
HH troscopic, stars: early-type, stars: winds, outflows 



O 1. Introduction 

H . 

During recent years many OB stars have been discovered to be binary systems. Non- 
, thermal radio emission is observed to be produced by some of these binary stars. The 
non-thermal emissivity is thought to be due to synchrotron emission radiated by rel- 
. ativistic electrons. The electrons are accelerated up to high energies by strong shocks 
J> ' produced by the collision between the two radiatively driven stellar winds (Eichler & 
, Usov 1993). Several parameters of the system can be constrained by the synchrotron 
- emission, among them the mass loss rates from the primary and the secondary. Investi- 
gating the synchrotron radiation is thus necessary. We model the non-thermal emission 
for two colliding wind systems, Cyg 0B2 No. 8A and Cyg 0B2 No. 9, and compare the 
obtained results with the observations. 
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2. Modelling and comparison with the data 

The two colliding winds are separated by a contact discontinuity. Its position (which, 
in our model, we assume to be coincident with the two shocks) is defined as in Antokhin 



^ et al. (2004). The electrons are accelerated at the shock. We follow them as they advect 
" " ' away and cool down due to adiabatic and inverse Compton losses along the post-shock 
streamlines. The momenta follow a modified power-law distribution. The synchrotron 
emissivity from the relativistic electrons is calculated along the post-shock streamlines 
in the orbital plane. The Razin effect is included. The third dimension is recovered by 
rotating the orbital plane along the line which connects the two stars. We also include 
free-free emission and then calculate the fluxes and spectral indices at different orbital 
phases using Adams method (jAdam 1990p . For Cyg 0B2 No. 8A the parameters are 
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Figure 1. Flux in the radio band at 6 cm for Cyg OB2 No. 8A: on the left the observations 
from Blomme et al. (2010), on the right our simulated results. Periastron is at phase ~ 

provided by De Becker et al. (2006), for Cyg 0B2 No. 9 the orbital parameters are 
provided by Naze et al. (2010), stellar parameters by Martins et al. (2005) and wind 
parameters by Vink et al. (2001). 

Cyg OB2 No. 8A results. For Cyg 0B2 No. 8A (see lBlomme et al. 2010|) the radio 
data at 3.6 and 6 cm are obtained with VLA, the X-ray data with XMM and ROSAT. 
Variability that is locked with the orbital phase is observed in both radio and X-rays. The 
X-ray and the radio light curves are anti-correlated due to different formation regions. 
The radio formation region is far out in the wind, along the contact discontinuity, while 
the X-rays are formed much closer to the apex of the contact discontinuity. The model 
predicts phase-locked radio variability which is consistent with the observations (see 
Fig. [T|), even if the phases of the flux maximum and minimum do not agree. 

Cyg OB2 No. 9 results. IVan Loo et al. 2008l studied the observed VLA radio fluxes 
of Cyg 0B2 No. 9 at 3.6, 6, and 20 cm. They found a 2.35 yr period from the data. 
A preliminary 6 cm light curve from our modelling shows variability in the radio flux 
linked to the orbital period that is the fingerprint of non-thermal radiation. Compared 
to the observations, the theoretical fluxes are much too high and the maximum occurs 
too early. 



3. Future work 

To improve the current results for Cyg 0B2 No. 8A and 9, we need to include the 
orbital motion and solve the hydrodynamics equations. A more detailed study is also 
necessary to better determine the star and wind parameters of the binary components 
and to investigate the porosity /clumping problem. 
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